Variational speed selection for the interface propagation in superconductors 
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We study the interface propagation in superconductors by means of a variational method. We 
compute the lower and upper bounds for which the planar front speed propagation is valid. To take 
into account delay or memory effects in the front propagation, an hyperbolic differential equation is 
introduced as an extension of the model. 
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In recently years the magnetic field penetration or 
its expulsion from Superconducting samples have at- 
tracted the attention of different research groups. The 
understanding of the different particularities of this phe- 
nomenon has been a major challenge. 

In RefU, the authors have obtained the interface speed 
by using heuristic methods such as Marginal stability hy- 
pothesis(MSH) and Reduction order. 

In this paper, we aboard the interface propagation 
speed from the variational point of view. The variational 
speed selection or BD method was proposed by Benguria 
and Depassieii^i^i^ in order to study the reaction-diffusion 
equations. Using a trail function g{x) in the procedure, 
one may find accurate lower and upper bounds for the 
speed c. The function g{x) must satisfy that g{x) > 
and g'{x) < in (0, 1). Only if the lower and the upper 
bounds coincide can the value of c be determined without 
any uncertainty. 

Our start point are the Ginzburg-Landau equations^, 
which comprise a coupled equations for the density of 
superconducting electrons and the local magnetic field. 
In order to describe the evolution of the system between 
two homogeneous steady states, we assume a SC sample 
embedded in a stationary applied magnetic field equal 
to the critical H — He. The magnetic field is rapidly 
removed, so the unstable normal-superconducting planar 
interface propagates toward the normal phase so as to 
expel any trapped magnetic flux, leaving the sample in 
Meissner state. Also, we have considered that the inter- 
face remains planar during all the process. 

To take into account the delay effect in the interface 
propagation, due to, for example, imperfections and non- 
homogeneous superconducting properties in the material, 
we have included the delay time t and indeed introduce 
the hyperbolic differential(IID) equation. This type of 
equation has been recently applied in biophysics to model 
the spread of humans^, bistable systems^, forest fires^ 
and in population dynamics^. 

Traveling wave solutions. We are interested in find- 
ing traveling wave solutions for our model. To start 
we use the one-dimensional time-dependent Ginzburg- 
Landau equations (TDGL), which in dimensionless units^ 



dtf ^- dlf-q'f + f~ /^ 
adtq - dlq - fq, 



(1) 



where / is the magnitude of the superconducting order 
parameter, q is the gauge-invariant vector potential (such 
that h = dtq is the magnetic field), a is the dimension- 
less normal state conductivity (the ratio of the order pa- 
rameter diffusion constant to the magnetic field diffu- 
sion constant) and k is the Ginzburg-Landau parameter 
which determines the type of superconducting material; 
K < 1/V2 describes what are known as type-I supercon- 
ductors, while K > I/V2 describes what are known as 
type-II superconductors. 

In our analysis we will search for steady traveling waves 
solutions of the TDGL equations of the form f{x,t) = 
s{x — ct) and q{x, t) = n{x — ct), where z = x — ct with 
c > 0. Then the equations become 



Szz + CS2 — n s + s 



0, 



acuz — s n — 0, 
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(2) 



Vector potential q = 0. In this section, we assume 
q = for the TDGL equations, 



dtf ^\dlf + .f- f. 



(3) 



Then, there exists a front / = s{x — ct) joining / = 1, 
the state corresponding to the whole superconducting 
phase to / = the state corresponding to the normal 
phase. Both states may be connected by a traveling front 
with speed c. The front satisfies the boundary conditions 
lims^-^oo / = 1, lim^^cx; f = 0. Then Eq.Q can be writ- 
ten as. 



Szz + CK^ Sz +dk{s) =0, 

where is given by 5"^ = s{l — s"^) . 



(4) 
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c > 2 



./^ /o [n\l-nf{l-n){2-2n)]'''ds 
j',{\-nfds 



Finally, for the Eg. pTj) we arrive to c > 1, which is a 
better lower bound than the \/2 predicted by the MSH 
method. 



II. FRONT FLUX EXPULSION WITH DELAY 



FIG. 1: Predictions for the speed. The continuous curve is 
the lower bound predicted by the BD method given in Eq.((8l. 
The dashed curve is for MSH. 



We define p{s) = —ds/dz, where the minus sign is 
included so that p is positive. One finds that the front is 
solution of 



ds 



(5) 



with p (0) = 0, p (1) = 0, p > in (0, 1). 

Let g be any positive function in (0, 1) such that h = 
—dg/ds > 0. Multiplying Eq.fO by g{s) and integrating 
by parts between s = and s = 1 and taking into account 
hp + (g-Sk/p) ^ 2 y/ g h , we obtain that, 



c > 



^ 1 
{gh^)i ds/ / g ds, 



(6) 



As a trial function we have chosen g{s) = (1 — s)'^. Then 
one finds that, 



c > 



'j^sil~s)Hl-s')i2-2s) 



1/2 



ds 



J,{l-s)^ds 



after integration the speed is given by, 
3 



c > 



64 /fc 



124 + 37V2 log(3-2V2 



(7) 



(8) 



Notice that c < 2/k, where 2/k is the result obtained by 
using the MSH method. In Fig.l, the graphic shows that 
for values k > 1.4 the MSH speed tends to the BD value, 
but for K < 1.4 the variational speed selection provides a 
better lower bound. 

Vector potential q = 1 — f. For a set of parameters^ 
K = I/V2 and a = 1/2, we have that s{z) + n{z) = 1, 
then Eq.Q takes the form. 



(9) 



With this in mind, we look for solutions of the form 
s{z) = 1 — n{z). Proceeding as in Eq.® we have that. 



c> 2V2 [ (ghd)^ ds/ [ g ds, 
Jo Jo 



(10) 



An import feature phenomena is the existence of a de- 
lay time. In systems with interface propagation, this can 
be taken into account by resorting to the hyperbolic dif- 
ferential equation seen in Section I, which generalizes the 
parabolic equation. The aim of this section is to study 
the interface speed problem in superconducting samples 
by means of the HD equations. 

Our starting point is the HD equation. 



9^ u 



d u 
~dt 



d"^ u 
9^ 



fin) 



dfju) 
dt ■ 



(12) 



In the absence of a delay time (t = 0), this reduces to 
the classical equation ut = Uxx + f{u)- 

Vector potential q — 0. Taking into account the Eqs.((T]) 
and (|12p we can write the following expression. 



(13) 



where = s(l — s^). 

It has been proved^i^i^, that Eq. has traveling wave 
fronts with profile s{x — ct) and moving with speed c > 0. 
Then we can write Eq. p^ as follows, 

{I - ac^) s,, + c[n'' - a:S' {s)]s, + :Sk{s) = 0, (14) 

where z = x—ct, a ^ t, ~ ^, and with boundary 



conditions limz 



0, liniz^-c 



1. and < in 



(0, 1); Sz vanishes for z ±00. 

For the variational analysis we define p (s) = —Sz with 
p{0) = p (1) = and p > in (0, 1). Then the Eq.fTil) 
may be written as 

(1 - a c^) p ^ - c _ a y (s)] p + ^^(g) = 0. (15) 
ds 

Multiplying Ea. (|15p by g/p where g is an arbitrary posi- 
tive function and integrating by parts, we have that 



/\[i_4y']ds= f\{l-ac')hp+^]d 
Jo ' K Jo P 



s(16) 



where we have used the relation 



(1 - ac^)hp + ^ > 2V1 - ac^ y/ghdk , (17) 



P 



3 



and h = — g' > 0. 



a c 



Jo 9{>^'^ -a^')ds 



The maximum is attained for a g. Thus, the expression 
for the velocity is given by 



c> 2k.- 



[/2 + A ally- 



(19) 



h = 



g{K^ ~ aS')ds, (20) 



Notice that if the delay time is neglected a = 0, this 
reduces to Eq.®. 

The lower bound. To compute the lower bound we start 
with the trial function given by g{s) = (1 — s)^ and the 
expression for which both are substituted in Eq. (|19p . 
Then, 



(21) 



h = / [2n{n^ - l){n-lf] ^ dn, 
Jo 

( {^~nf [k^ ~ a[l~in^)\ dn. 
Jo 

from Ea. lfTU)) we have that, 
c > 2k — 



J 



where 



J 



[1 + 4K2a J2]2 

15 [124 + 37%/2 log (3 - 2^2)] 
64(10F - 7 a) 



(22) 



(23) 



The upper bound. The upper bound can be computed 
by using the Jensen's inequality^, 



< 



Jq fi{s) a{s) ds 



(24) 



where /i(s) > and a{s) > 0. We define ^{s) = g{K^ - 
a^') and a(s) = ^h/giK^ — a^')^. Then we can write 



J^\gh^y/'ds 
/o^5(«^ - ad')ds 
where 



< 



[hS/{K^-ar)] ds 
Jo gi^^"^ - aS')ds 



1/2 



(25) 



n\l-n^i2-2n)] , 
-ds — I -r^ ^ an. 



lo (K^-a^') Jo [k^ ~ a{l - 3n^)] 

/ g{K^-a^')ds= {l-nf[K^-a{l-3n'^)]dnl26) 
Jo Jo 

then we have that, 
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FIG. 2: Time-delayed interface propagation speed for q = 
0. The plot of the lower and upper bounds from variational 
method as well as the speed proposed by linear stability are 
shown. 



c < 2k- 



I 



[1 + 4K2aJ2]2 



where, 



and, 



a^ 



X= (1/3) ih+h)' 



(27) 



(28) 



(2a + k"^) log 



2a + K^ 



- (5a + 2k2 



2(k2 - a) 5 (2a + k^) arctan ( , /-^J^ 



V3a5/2 



(29) 



The Eq. ((27|) gives a better upper bound than the one 
predicted by linear stability- i.e, c < Cma-x = l/\/a- 

In Fig. 2 we have plotted the results of the BD method 
given by Eqs. iP^ and ([77|) as well as the bound proposed 
by linear stability(LS) methodology. The interface speed 
propagation can be predicted in a precisely way by using 
our trial function. On the other hand the difference with 
linear stability result is notable. 

Vector potential q = 1 — /. Taking into account the 
Eqs.lH) and IT^ we can write the following expression. 



T d^f 19/_5V 1 TdS_ 
2 9*2 + 2 at - dx^ ^ 2^^ 2 dt' 

where 5^ = s2(l — s). 

Then we can write Ea. ((30)) as follows. 



(30) 



{l-ac^)s,, + c[K^-ad'{s)]s,+dk{s)=0, (31) 

where we have assumed = (1/2)J? and a — t/2. 
The expression for the velocity is given by 



VI - ac2 



> 



g.^/o Ioi9hdV/'ds 

J^g{l-2a^')ds 



(32) 
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FIG. 3: Time-delayed interface propagation for q = 1 — /. The 
plot of the lower and upper bounds from variational method 
as well as the speed proposed by linear stability are shown. 



Proceeding as in Eg. p^ . we get the following expres- 
sion, 



> 2V2 



{3i + 8a3f) 



^2^l/2 



(33) 



where 



li = [ ^/gh^ds, 32 = ( g{l-2a^')ds, (34) 

The lower bound. As mentioned before, one may obtain 
lower bound for the interface speed by means of our trial 
function g{n). Taking into account Eqs. ([M)) . the integral 
functions can be written as, 



31= [n'^{2-2n){l-nfY dn, 
Jo 

{1-nf [l~2a{2n- 371^)] dn. (35) 

Then the velocity takes the form, 

c > 75 (15 -a)- 1/2, (36) 

The upper bound. To compute the upper bound we 
have used the expression Eg. ((25)) but with ^ = s^(l — s) 



and = 1/2, then 



< 



hd 



J^g{l-2ad')ds 



(37) 



/ (1 



2 a;?') "^'^ L [n2(2-2n)(l-n)]' dn, 



1 .1 
g{l-2ad')ds = / (1 - n)^ [l - 2 a (2 n - 3 n^)] dr(38) 
^0 



After integrating and do some algebra, the expression 
for the velocity is given by 



> 2V2 



03 



1/2 ' 



(39) 



where. 



1/2 



„ 01 2 arctan a arctan 4a „ , , , 

= — ^7= , + — J—- +02] , 40 



a= ^al^--aj/i3^a), (41) 

/?! = 80^ + 6a -I- 9, 
p2 = 2V^ [-3(3 + 4a) + (3 -f 2a) log(l + 2a)] . 

In Fig. 3 we have plotted the results of the BD method 
given by Eqs. ((36|) and p9)) as well as the bound proposed 
by LS method. The interface speed propagation can be 
predicted in a precisely way by using this trial function. 
On the other hand the difference with LS result is no- 
table. 

Conclusion. Throughout this work, we have performed 
analytical analyses on the superconducting-normal inter- 
face propagation speed problem in parabolic and hyper- 
bolic equations. We have made use of the variational 
analysis to obtain the lower and upper bounds for the 
speed in each case. 
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